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Abstract: A detailed numerical analysis is performed to obtain the Hawking spectrum 
for charged, massive brane scalars and fermions on the approximate background of a brane 
charged rotating higher-dimensional black hole constructed in [1]. We formulate the prob- 
lem in terms of a "spinor-like" first order system of differential wave equations not only for 
fermions, but for scalars as well and integrate it numerically. Flux spectra are presented 
for non-zero mass, charge and rotation, confirming and extending previous results based 
on analytic approximations. In particular we describe an inverted charge splitting at low 
energies, which is not present in four or five dimensions and increases with the number of 
extra dimensions. This provides another signature of the evaporation of higher-dimensional 
black holes in TeV scale gravity scenarios. 
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1. Introduction 

In a recent paper [1] an approximate background to model the gravitational field of a higher 
dimensional rotating black hole with a brane confined abelian charge was constructed. This 
background was coupled to massive, charged brane scalar and fermionic perturbations and 
the corresponding wave equations were separated to allow for the study of the associated 
Hawking radiation [2]. 

The detailed study of black hole perturbations is important in many different contexts 
such as astrophysics [3-7], cosmology [8] or from the purely quantum field theoretical 
point of view [9-11]. The particular construction in [1] was motivated by TeV gravity 
scenarios which contain extra dimensions [12-17], and the Standard Model fields confined 
to a 4-dimensional brane [18-20]. In such scenarios black holes may form [21-24] in the 
high energy collision of charged brane degrees of freedom so the brane charge of the black 
hole becomes important. In particular, recently, it has been shown numerically in four 
dimensions that black holes indeed form in ultra-relativistic collisions between Schwarzchild 
black holes [25-27] or between solitons [28]. These results support the early hoop conjecture 
type arguments. 
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In this paper we complete the study in [1] by performing a full numerical analysis 
which is not constrained by the approximations of small energy, mass, charge and rotation 
used in analytic approaches. This allows us to obtain exact numerical results in the full 
range of energies which is important even if some of the parameters such as charge, mass 
or rotation are small. We confirm all the low energy features found in the analytic study, 
in particular the weakness of discharge for typical QED like charges. We also present plots 
for larger masses which are relevant for typical TeV gravity scenarios at the LHC. 

The structure of the paper is the following: In Sec. 2 we briefly review the background 
fields and the wave equations for the scalar and fermionic perturbations with mass and 
charge. We formulate the scalar equation in a new "spinor-like" form which is particularly 
convenient for numerical integration. In Sect. 3 we provide a series expansions for the fields 
near the horizon, an asymptotic expansion at infinity and re-formulate the first order system 
in a more convenient form to extract the transmission factor. In Sect. 4 we present plots for 
transmission factors and the number flux spectra, confirm the low energy results, extend 
them to the full energy range and point out the main features with and without rotation. 
Finally in Sec. 5 we summarize the results emphasising the importance of implementing the 
new effects in black hole event generators [22,29-32] to perform phenomenological studies 
of TeV gravity with black hole production. 



2. The fields 

We want to study a charged rotating black hole background where the abelian (i.e. Maxwell) 
field is confined to a four dimensional brane. As derived in [1] an approximate effective 
background in 4 + n dimensions is characterized by the metric 



where 

A = r 2 + a 2 + Q 2 --^ T , S = r 2 + a 2 cos 2 B , (2.2) 

p, = 1 + a 2 + Q 2 , if we take horizon radius units (r# = 1 - see Sect. 4.1 of [1]); and the 
Maxwell field is 

A a dx a = -Q^ (dt - asm 2 0#) . (2.3) 

Here Q is the electric charge of the black hole and a its oblateness parameter which is 
directly related to the amount of angular momentum [33]. If now we consider coupling 
other quantum fields of various spins, Hawking radiation is present [2,34-38] and the fields 
are thermally emitted from the hole which evaporates progressively [39-49]. The various 
fluxes of particle number N, energy E, angular momentum J and charge Q, are given by 

d 2 {N,E,J,Q} _ 1 ^ J {l,w,m,q} T (4+n), m ,„ A) 
^ " 27 ^ ^ exp(fi/31r)±l Tfc K °' 9 ' Q) ' {2A) 

j=\s\m=-j 
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where k = {j, m} are the angular momentum quantum numbers of the partial wave; uj, fi, q 
are the energy, mass and charge of the particle respectively, Co = co — (ma + qQ)/(l + a 2 ) 
and 

(n + l) + (n-l)(a 2 + Q 2 ) 
TH = 4vr(l + a 2 )r H (2 " 5) 

is the Hawking temperature. The term containing the exponential in (2.4) is the so called 
Planckian factor. The transmission factor T^ 4+n ^ is the fraction of a wave incident from 
infinity which is transmitted down the horizon and is purely ingoing at the horizon. This 
factor is obtained by solving the wave equations for the particular field with such ingoing 
boundary conditions at the horizon. For the background (2.1), (2.3) separation of variables 
yields the following radial and angular equations [1]: 



Massive charged scalars: 



where 



K = oj(r z + a 2 ) - am - qQr (2.7) 
The boundary condition at the horizon is [1,50] 



U = fi 2 r 2 + A CJ - jm + oj 2 a 2 — 2auim . (2i 



R = x {K So (1 + ...) (2.9) 

with x = r — 1, K± = u>(l + a 2 ) — am — qQ, and do = n + 1 + (n - 1)(1 + a 2 + Q 2 ) 
is the leading order coefficient of the expansion of A in powers of x. The angular 
eigenvalue A C) j jTn is determined from the angular equation 

|( sm ^) + ( c W*--^ + A c ,„„)s = (2.10, 



sin 9 d6 

with c 2 = a 2 (uj 2 — fJ 2 ), by imposing regularity of the solution at cos 9 = ±1. For 
a = we have the closed form Aoj, m = j(j + 1). 

Eq. (2.6), can be written as a first order system of differential equations. This will 
be useful to perform the numerical integration using a method similar to that for 
fermions. Since there is no unique way of reducing the second order equation to a 
first order system, we take advantage of the extra freedom to construct a spinor-like 
object with a conserved Wronskian and, simultaneously, an asymptotic behaviour at 
infinity which gives the transmission factor straightforwardly. It is then possible to 
show that a convenient choice is 

P ± » = #(**<), (2.1D 
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where in principle k can be an arbitrary constant but we set it to the momentum of 
the partial wave k = \j uS 1 — fi 2 . So the second order equation (2.6) is replaced by 
the first order coupled system 1 

d~P 

s =M,(r)P s (2.12) 



dr 



where 



^M^-i + + (2 ' 13) 

and &i are the Pauli matrices. Now, using (2.12), conservation of the Wronskian is 
easily checked: 

T r { P ^s)=l(\P +lsl \ 2 -\P- lsl \ 2 )=0. (2.15) 

The choice k = y 'uS 1 — n 2 ensures that P±q picks respectively the outgoing/incoming 
part of the wave at infinity (see Sect. 3.2). 

Massive charged fermions: For fermions, the radial equation obtained in [1] is already 
in the form (2.12) with 

, «■ , s A „ tir „ K ^ , . 

Mi r = — a x - -^a 2 + i-ra 3 . 2.16 

2 ^2 A2 A 

So P 1( / 2 obeys the same Wronskian relation (2.15) as does (2.12). Again, the incoming 
solution at the horizon takes the form 

Pi ~x _i ^f(ao + ...) (2.17) 

2 

with ao a constant spinor. The angular eigenvalue is obtained from the system of 
angular equations 



d „ / . n m \ 1 

— + 2s aw sin H — cot ( 

d6 V sinOJ 2 



S-s = {2s\ + an cos 0)S S , (2.18) 



where s = ±1/2. Once again, in general, the eigenvalues are obtained by imposing 
regularity of the solution at cos# = ±1. When a = 0, A = j + 1/2 with j a positive 
semi- integer. 

Electromagnetic and other perturbations: Regarding electromagnetic perturbations, 
it is known in four dimensions that they couple to gravitational perturbations for the 
Kerr-Newman black hole (see for example Sect. Ill, Chapter 11 of [51] and refer- 
ences therein). Similarly, in the higher dimensional case, we would expect them to 
couple to gravitational modes on the brane. However, in the limit of small charge the 



1 Here s is the spin which we leave arbitrary since the same type of equation will hold for fermions. 
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perturbations should approximately decouple. This is indeed the case and an approx- 
imation scheme was developed by Dudley and Finley [52]. It amounts to considering 
separately one perturbation (either electromagnetic or gravitational) while setting 
the other to zero on the fixed background. This approximation was used for example 
in [53] and [54] to compute quasinormal modes. In [54] this was compared to other 
methods to confirm the validity of the approximation for small Q (the special case 
J = was used) . The approximate second order wave equation for a perturbation of 
spin s is [54] 



a 1-s d 
A s — 

dr 



dr 



+ 



' dA ( dK 

K — is——K + A I lis- A 

dr \ dr 



R = 0, (2.19) 



where K is the same as in (2.7), but does not contain the particle charge g-term. 
This correctly reduces to the exact result for scalars and fermions when a back- 
ground charge is present (if K contains the g-term in (2.7)) and it describes the 
electromagnetic or gravitational perturbations approximately, for small Q. 

An important feature of electromagnetic and gravitational perturbations, compared 
to scalars and fermions, is that they are electrically neutral. No electric coupling 
means that qualitatively, not much will change compared to the case of no back- 
ground charge. Specially for small charges we see from (2.19) that the charge of the 
background only enters through the Q 2 term in A. This affects mostly the Hawking 
temperature which in the small Q limit will simply rescale the flux curves without 
much difference in shape. In fact as noted in Fig. 2 of [1] for scalars and fermions, 
the effect of a small background charge on neutral particles is indeed small on the 
transmission factors and the flux curves are simply rescaled by the different Hawk- 
ing temperature in the thermal factor. So qualitatively nothing changes for neutral 
scalars and fermions which obey (2.19) so we would expect the same for higher spins. 
Therefore we will not present numerical results for the electromagnetic or gravita- 
tional field, since they reduce to well studied cases (see e.g. [40,44,47]) both qualita- 
tively and in terms of implementation (the following constant shift a 2 — > a 2 + Q 2 in 
A is sufficient). 

Furthermore, if we assume electroweak symmetry is not restored outside the black 
hole and that the electrically charged weak vector boson W and the neutral Z pro- 
vide a good effective description of the weak degrees of freedom, it is tempting to 
guess that (2.19) holds similarly for those perturbations (with K containing the elec- 
tric g-coupling). This is because the black hole background can only be electrically 
charged (or colour charged) so the backrgound values of the weak field perturba- 
tions vanish. Then we would expect (2.19) to be exact since there is no reason for 
the weak field perturbations to couple to the linearized gravitational perturbations. 
This is in contrast with the equations for electromagnetic perturbations where terms 
linear in the gravitational perturbations arise from linearising bilinears in the grav- 
itational/electromagnetic fields around their background values. For weak W and 
Z field perturbations (as for scalars and fermions) such gravitational terms can not 
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be present because even if they exist before linearisation, when evaluated on the 
background for the W and Z fields they are identically zero. 

Furthermore, because the W and Z fields are massive, they are described by a complex 
or a real Proca field respectively, which is an extra complication. 

Alternatively, if electroweak symmetry is restored in the region outside the black 
hole 2 , then we have to use the fundamental weak gauge fields associated with the 
SU(2)l x U(1)y sector of the Standard Model (instead of the electromagnetic, the 
W and the Z fields). 

Due to these extra complications, the detailed study of other vector perturbations 
will be treated elsewhere. 



3. Numerical methods 

In this section we present the methods used to reduce the linear systems of equations at 
hand to initial value problems which are more convenient for numerical integration. 

3.1 Near horizon expansions 

The boundary condition at the horizon is most easily implemented through a series ex- 
pansion. This allows for a high precision initialisation of the radial functions slightly away 
from the horizon to avoid numerical difficulties associated with the coordinate singularity. 
The expansions we need are 



R 



Pi 



£ Ct r- 
m=0 

+00 

a ™ (v^)' 



X > 

m=0 

+00 

x a 

m=0 



(3.1) 



Note that R can be used to initialise Po- By inserting into the wave equations (2.6) 
and (2.12) respectively we obtain the following recurrence relations 



a 



So 

a = 1, a m = 



a 



m— 1 




m(m + 2o)(5q 
a m = (N - 5 (m + 2a))- 1 



(m + a)J 7m + ^ k ( k + ^m-k + a k a m - k ) 



k=0 



m—1 

3=0 



m > 1 



(3.2) 



where a choice of normalisation was made, when setting ao and ao. The various coefficients 
are defined in appendices A.l and A. 2. Using expansions (3.1) we have initialised P s at 



2 This should be the case if the black hole size is smaller than the electroweak breaking scale which is 
typically 1/mw, the inverse mass of the W. 
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x = 0.1 by truncating the series at eighteenth order. A first estimate of the numerical error 
can be made by modifying this choice (we have used x = 0.05 and x = 0.01 as a check). 3 



3.2 Far field expansions 

Once the radial function is initialised, numerical integration routines can be used to prop- 
agate the solution away from the horizon according to (2.12). When sufficiently away from 
the horizon, the transmission factor can be evaluated by comparing the numerically propa- 
gated solution with its asymptotic form at large r. An asymptotic expansion can be found 
in the form 



e qr r -l 



m=0 



s r -m 



if we expand 



m=0 

and equate (2.12) order by order. The leading behaviour is 

p = y a (°ut) e nyvd+ + Y^e^y-^dJ , 
where y = kr, Yj out ^ and Y^ are constants, 



e k~ a k 



e = ~qQ + w(l + a 2 + Q 2 )5 nfi a = |(1 + a 2 + Q 2 )5 nfi , 



and 











u)+k 



1 



(3-3) 
(3.4) 

(3.5) 

(3.6) 
(3.7) 

(3.8) 



We can now factor out the dependence at infinity so that the leading asymptotic form 
for the upper (lower) component of the spinor becomes Y} out ^ (Yj m ^ ) respectively. This is 
achieved by performing a rotation on the spinor P s such that it eliminates a fixed number 
of subleading terms in the asymptotic expansion (3.4) (in practise we have eliminated the 
first two subleading terms). Then the new spinor Q s is related to P s through Q s = R S P S 
and the the system to integrate becomes 



dQ s 

dy 



(3.9) 



The explicit forms for A s and R s are given in appendix B. 

Finally, the transmission factor is computed from the definition by taking the limit 



T, 



lim 

r— >+oo 



1 - 



Q 



+s 



Qs 



Y 



(out) 



Y 



(in) 



(3.10) 



3 Throughout we have required an error s < 10 4 , for the transmission factors. 
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(±s for upper/lower component respectively) and an estimate of the error is obtained by 
varying the large r used in the limit. Furthermore, with the normalisation chosen in (3.2) 
we can evaluate the Wronskian (2.15) at the horizon and use its conservation to obtain a 
second expression 

kW s kW s 

T s = lim — * = — , (3.11) 



where 



W = K* Wi= . (3.12) 



By comparing the results from (3.10) and (3.11), we obtain another estimate of the nu- 
merical errors. Eq. (3.11) is particularly useful since it contains explicitly the zeros of the 
transmission factor in the numerator. 

To integrate (3.9), a code was written in C++ using the Gnu Standard Library (GSL) 
numerical integration routines. This was checked against an independent code in Maplell. 

3.3 Angular eigenvalues and angular functions 

To determine the transmission factors when the rotation parameter a is non-zero, it is 
necessary to solve the angular equations (2.10) and (2.18) numerically (no closed form is 
known for the angular eigenvalue when 0). 

Using (2.10) and (2.18), it is easy to show that all cases, except for massive fermions 
in a rotating background, are described by the following second order equation: 

Id/ „dS s \ f 99 „ „ (m + scosfl) 2 , \„ , . 

where for fermions A c j jTn = X 2 — a 2 uj 2 + 2aumi — \s\ — s. This equation, which describes 
spheroidal harmonics, has been studied extensively in the literature. Whenever we evaluate 
the result with a / 0, we adopt the method in Appendix D of [32] to obtain the angular 
eigenvalues. 

For massive fermions on a rotating background, Eq. (3.13) will contain an extra term 
linear in afidS s /d9 (see Eq. (4.29) of [1]) so the method of [32] cannot be applied. Nev- 
ertheless, since we are mostly interested in studying mass and charge effects, Eq. (3.13) 
allows us to obtain a representative set of cases 4 . 



4. Numerical results 



In this section, samples of numerical data of transmission factors were generated using the 
method presented in Sect. 3. From such data all interesting fluxes and distributions can 
be computed quickly. Most of the samples were generated up to w = 10, but some up to 
co = 5 to save computing time. We show plots with oj < 5 since the curves are very quickly 
stabilised for large oj (either to a constant or a suppressed tail). 

4 Note that except for some final plots (which are present for illustration purposes), all the results in 
Sect. 4 will have a = 0. 
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Figure 1: Scalar (left) and fermion (right) transmission factors for n = 2: The top plots are for 
variable [i and the bottom plots are for variable q. The first three partial waves are presented (note 
that a — 0, so waves with different m for the same j are degenerate). 

In sections 4.1 and 4.2 we focus mostly on results with the rotation parameter off 
and describe the main features for different charges and masses in the full energy range 5 . 
For illustration purposes, in section 4.2 we present some curves with typical rotation and 
typical charges which may be relevant for TeV gravity scenarios. 

4.1 Transmission factors 

In this section we display transmission factors for individual modes, confirm the earlier 
results obtained in [1] and extend them to the full energy range. We describe the main 
features of the plots which are relevant to the discussion of the fluxes in the next sections. 
We have checked that the approximate results based on the analytic approximations of [1] 
reproduce well the exact numerical results obtained with our method even at intermediate 
energies as claimed there. 

In Fig. 1 we present plots for the n = 2 case and a range of charges and masses, with 
id G [0,2]. 6 The top plots show the first three partial waves for scalars and fermions and 
[i = 0, 0.5 and 1. For scalars we confirm the strong suppression at the mass threshold for 
the j = partial waves, and the shift and suppression for higher partial waves [1]. For 
fermions the behaviour is similar, except that the first partial waves are not so sharply 
suppressed at threshold. The bottom plots show the same partial waves when the mass is 
set to zero, the background charge is set to Q = 0.4 and the particle charge varies between 

J The effect of rotation was studied previously in [39,41-44,47-49], and massive fermions without rotation 
and at low energies in [55]. Here we are mainly interested in charge and particle mass. 
6 The transmission factors asymptote to unity quickly, so this is the interesting region. 
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Figure 2: Scalar (left) and fermion (right) transmission factors for variable n: The top plots are 
for variable /i and the bottom plots are for variable q. Within partial waves with the same j, the 
curves are ordered from n = 1 to n = 6 from top to bottom (e.g. j = 1 curves indicated in the 
plots). 

q = 1 and q = — 1. Again we confirm, for both scalars and fermions, that negative charges 
are favoured in the full range of energies (all curves split following the same pattern as 
indicated for the first mode). This is because the transmission factor is the fraction of a 
wave incident from infinity that is transmitted through the horizon, so we would expect 
the Coulomb attraction to favour such negative charges as observed. This feature will be 
important to understand the behaviour of the fluxes, in the next section. 

Figure 2 shows the variation with n. The top plots are for \x = 0.5 and the bottom 
plots for q = 1 and Q = 0.4. The general tendency is for the transmission factor to be 
suppressed with n. The exception is in the low energy region, when the charge is non-zero, 
where the tendency is inverted for fermions, whereas for scalars in the superradiant region 
the variation with n is small (this agrees with the results in [1]). 

4.2 Fluxes 

In this section we present plots for the particle number flux in the full energy range, summed 
over partial waves (we have included the first ten j-partial waves). This quantity is sufficient 
to illustrate the new effects, since for example the power flux curves are qualitatively similar 
and the charge fluxes are proportional to the corresponding number flux. Furthermore, the 
angular momentum flux (when rotation is present) has been studied in detail before [39,41- 
44,47-49] and we will observe that even when the new effects are present, the contribution 
from rotation affects the spectrum in a similar fashion as in those earlier studies. 

In Fig. 3 we present examples of the number flux for non-zero mass, when the charge 
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Figure 3: Scalar (left) and fermion (right) number fluxes for n = 4 and zero charges: The top 
plots show fi — 0.5 and the contributions from each partial wave to the total flux. The bottom 
plots show variable [i and variable n = 5, . . . , 2. For each /i the curves are naturally order in n from 
top to bottom, n = 5 and n = 4 are indicated for // = 0. 

and the rotation are set to zero. The top plots show n = 4 and \x = 0.5. Note that typical 
values of fj,, for Standard Model heavy particles such as the top quark in TeV gravity 
scenarios, 7 range from 0.1 to 0.5. We also indicate the contributions from the first few 
j values to the total flux curve. Similarly to the transmission factors, the main feature 
is a sharp suppression at threshold. The area under the curves is larger for scalars than 
fermions, which agrees with earlier studies (see for example [43,48]). The bottom plots 
show three values of the mass and various n values. We confirm the conclusion of [1] that 
the area under the curves is suppressed as \i increases and the suppression at threshold is 
smooth both for scalars and fermions. The error from using the \x = curve with a sharp 
cut at the mass is therefore large (most notably for fermions). Regarding variation with 
n, it is opposite to the tendency for the transmission factors so the n dependence of the 
Planckian factor dominates the magnitude. 

Figure 4 shows several cases of non-zero charges. We have kept Q = 0.6, which is a large 
value (see Sec. 3 of [1]) so that all effects can be seen easily. Similarly we show q in the range 
[—1,1]. We use n = 4 as a representative case. The top plots show the total flux for the two 
extreme cases q = 1 and q = — 1 together with the first few partial waves contributing. The 
first striking observation is the confirmation that for all partial waves there is a region at 
low energy where charging up is favoured (i.e. the curve corresponding to negative charge 
is higher) and then another (dominant) region where discharge is favoured (the curve with 

7 /i is in horizon radius units and 1/rsr is typically in the range 200 GeV — 1000 GeV for TeV gravity 
scenarios at the LHC. 
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Figure 4: Scalar (left) and fermion (right) number fluxes for n = 4, variable q, and Q = 0.6: 
The top plots show two opposite and large |g| = 1 cases to illustrate the charge splitting, together 
with the first three partial wave contributions. The bottom plots show the variation of the curves 
between these two large charges. 
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Figure 5: Scalar and fermion number flux asymmetries: Both plots show curves for the difference 
in number fluxes between positively charged and negatively charged particles for two values of \q\Q. 
The curves are naturally ordered in n (some cases are labelled) from n = (curve with the lowest 
maximum) to n = 6 (highest maximum). 

positive charge is higher). It is also clear that if we integrate over the curves discharge 
is always favoured as expected. The bottom plots show a similar behaviour for a range 
of intermediate charges. Another interesting point is that the splitting at low energies is 
larger for fermions than for scalars. 

The inverted splitting at low energies is a direct consequence of the extra dimensions. In 
Fig. 5 we show the difference in number flux of positively and negatively charged scalars and 
fermions when n = 0, . . . , 6. The left plot shows a typical QED coupling of \q\ = Q = 0.1 
and the right plot a QCD like coupling of \q\ = Q = 0.3. Note however that we are 
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Figure 6: Scalar and fermion asymmetries for a — 0.9: The left plot shows number fluxes for 
scalars and fermions with positive and negative \q\ — 0.3 for n — 4. The difference between positive 
and negative \q\ curves is also shown in the same plot. The right plot shows the difference for scalars 
and a range of n's. 

dealing with an abelian theory so the latter is only indicative of the magnitude of the effect 
for QCD. Prom this figure it is now clear that the splitting is controlled by an interplay 
between the transmission factor (which prefers negative charges) and the Planckian factor 
(which prefers positive charges). For n = and n = 1, the splitting is always positive 
so the Planckian factor dominates. However as n increases, the transmission factor starts 
dominating at low energies and for all n > 2 we have the observed inverted region (where 
the curves are negative). Another interesting feature of Fig. 5 is that the plots on the left 
have exactly the same shape as the ones on the right. This is not surprising if we note 
that for qQ small we can expand the fluxes perturbatively around qQ = and since \qQ\ 
is 0.01 and 0.09 respectively, we would expect the perturbation to be dominated by the 
linear term so the difference is proportional to \qQ\- 

Finally Fig. 6 shows some cases with a rotation parameter a = 0.9 (the typical order 
of magnitude for a TeV gravity scenario rotating black hole) and the typical QCD charges 
\q\ = Q = 0.3. The left plot shows the split flux curves for the QCD case both for scalars 
and fermions and the difference between the two. The right plot shows the difference 
curves for scalars and a range of n's. Qualitatively, the splitting of the curves when 
follows the same pattern as Fig. 5. The main differences are: the oscillations, which are 
due to the contribution of higher partial waves when a ^ (they are responsible for 
shedding the angular momentum of the black hole [39,41-44,47-49]); and the shift of the 
spectrum towards higher energies, which is again a well known effect of rotation related to 
the contribution of partial waves with larger j. It is interesting to note that the oscillations 
persist for large n in the right plot, which is not true for the flux plots where they tend to 
be smoother [39,41-44,47-49]. 

5. Conclusions 

In this paper we have performed a full numerical analysis of Hawking radiation for mas- 
sive, charged scalars and fermions on an effective higher dimensional rotating black hole 
background with charge. 
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In the first part we have re- formulated the problem in a convenient form and obtained: 
i) a series expansion to implement the boundary condition at the horizon, Eq. (3.1), and the 
asymptotic form at infinity, Eq. (3.3); and ii) a"spinor-like" first order system of differential 
equations both for scalars and fermions, Eq. (3.9). 

In the second part we have shown a selection of plots to illustrate each effect. The 
main results are: 

• We confirm and extend all the conclusions in [1] regarding massive particles to the 
full energy range, with all contributing partial waves. In particular we obtain the 
case of massive brane fermions for n > 0, which was not studied before in the full 
energy range. The main difference is that for fermions the suppression is not so sharp 
at the threshold energy. Since the typical mass parameters of heavy Standard Model 
particles for TeV black hole scenarios can go up to ~ 0.5, this is an important effect. 

• Regarding charges, we have confirmed the splitting of the fluxes between positive and 
negative charges in the full energy range and showed that discharge is always favoured. 
The most interesting feature we have found is the inverted charge splitting at low 
energies which, we have shown, is a new effect due to the extra dimensions for n > 2. 
So, even though electric discharge may be small in TeV gravity black hole events [1], 
this splitting will still be present and it may be possible to reconstruct it if such 
events occur in future experiments. For QCD charges, the splitting should be even 
larger but a non-abelian analysis will be necessary to determine which observables 
will display it. 

These conclusions remain qualitatively the same with rotation, which affects the 
fluxes in ways that were observed in earlier studies. 

To summarize, the methods we have described can be used in the full energy range to 
implement an improved model of the Hawking evaporation in black hole event generators 
with non-zero masses and charge asymmetries. 
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Appendices 

A. Expansion coefficients 
A.l Scalars 

The expansion coefficients we need are denned by 

+00 



A = x ^ S. 

m=0 

+00 

K 2 -AU=J2° m 



m=0 
m— 1 



7m 
Ira 



^2(k + a)a k 5 m _ k 

k=0 

(m + a)a m 5 + j m 



It can be shown then that 



5 = n + 1 + (n - 1) (a 2 + Q 2 ) 
i(n-l) (l + a 2 + Q 2 ) 



5 1 =1- 



<5 2 

^m+l 
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(n 2 -l) (l + a 2 + Q 2 ) 



-(1 + p m+ i)S v 



, m > 3 



where 



n - 2 



P2 
Pm+l 



1 



m + 2 



and 



where 
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(j 2 = 2K*w + (2w - gQ) 2 - U Q S 1 - U160 

(73 = 2u(2u - qQ) - U 5 2 - UxSx - U 2 5 
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a m = -Uo5 m -i - U\5 m -2 - U 2 6 m -3 ,m>5 



Uq = A + uj 2 a 2 — 2auim + p 2 



Ux = 2p 2 
U 2 =p 2 
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A. 2 Fermions 

Similarly to the scalar case define 



2AMi(r) = ^N m (V*T 

D 

+00 

m=0 
m— 1 

b 2 m = ^ 25 ™-i0' + «) a 2j 



A5 



j=0 
m—1 



b 2m +i = ^ S m -j(2j + 2a + l)a 2 j+i , 

j=0 



The matrices we need are 



N = 2iK*a 3 
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where <5j are obtained from the following expansion 
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(A.20) 
(A.21) 
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(A.24) 



(A.25) 



by fixing a certain order of truncation and expanding the square root in powers of x up to 
the given order. 

B. Matrices 

In the main text we have used the following matrices: 
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c m are coefficients such that the corresponding powers in the asymptotic expansion of (B.5) 
are cancelled; 
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and now 



A = y 2 + k 2 (a 2 + Q 2 )- 



(1 + a 2 + Q 2 ) A: n+1 



,,n-l 



References 

[1] M. O. P. Sampaio, Charge and mass effects on the evaporation of higher- dimensional 
rotating black holes, JHEP 10 (2009) 008, [arXiv: 0907. 5107]. 

[2] S. W. Hawking, Particle creation by black holes, Commun. Math. Phys. 43 (1975) 199-220. 

[3] T. Regge and J. A. Wheeler, Stability of a Schwarzschild singularity, Phys. Rev. 108 (1957) 
1063-1069. 

[4] S. A. Teukolsky, Perturbations of a rotating black hole. 1. Fundamental equations for 
gravitational electromagnetic and neutrino field perturbations, Astrophys. J. 185 (1973) 
635-647. 



-17- 



[5] W. H. Press and S. A. Teukolsky, Perturbations of a rotating black hole. II. Dynamical 
stability of the Kerr metric, Astrophys. J. 185 (1973) 649-674. 

[6] S. A. Teukolsky and W. H. Press, Perturbations of a rotating black hole. Ill - Interaction of 
the hole with gravitational and electromagnetic radiation, Astrophys. J. 193 (1974) 443-461. 

[7] S. R. Dolan, Scattering and absorption of gravitational plane waves by rotating black holes, 
Class. Quant. Grav. 25 (2008) 235002, [arXiv: 0801 . 3805]. 

[8] D. N. Page and S. W. Hawking, Gamma rays from primordial black holes, Astrophys. J. 206 
(1976) 1-7. 

[9] S. W. Hawking, Information loss in black holes, Phys. Rev. D72 (2005) 084013, 
[hep-th/0507171]. 

[10] C. Rovclli, Loop quantum gravity, Living Rev. Rel. 11 (2008) 5. 

[11] M. B. Green, J. H. Schwarz, and E. Wittcn, Superstring theory, vol. 1: Introduction, . 
Cambridge, Uk: Univ. Pr. ( 1987) 469 P. (Cambridge Monographs On Mathematical 
Physics). 

[12] I. Antoniadis, A possible new dimension at a few TeV, Phys. Lett. B246 (1990) 377-384. 

[13] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, The hierarchy problem and new 
dimensions at a millimeter, Phys. Lett. B429 (1998) 263-272, [hep-ph/9803315]. 

[14] I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, New dimensions at a 
millimeter to a Fermi and superstrings at a TeV, Phys. Lett. B436 (1998) 257-263, 
[hep-ph/9804398]. 

[15] N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phenomenology, astrophysics and 

cosmology of theories with sub-millimeter dimensions and TeV scale quantum gravity, Phys. 
Rev. D59 (1999) 086004, [hep-ph/9807344]. 

[16] L. Randall and R. Sundrum, An alternative to compactification, Phys. Rev. Lett. 83 (1999) 
4690-4693, [hep-th/9906064]. 

[17] L. Randall and R. Sundrum, A large mass hierarchy from a small extra dimension, Phys. 
Rev. Lett. 83 (1999) 3370-3373, [hep-ph/9905221]. 

[18] T. Flacke, D. Hooper, and J. March-Russell, Improved bounds on universal extra dimensions 
and consequences for LKP dark matter, Phys. Rev. D73 (2006) 095002, [hep-ph/0509352]. 

[19] N. Arkani-Hamed and M. Schmaltz, Hierarchies without symmetries from extra dimensions, 
Phys. Rev. D61 (2000) 033005, [hep-ph/9903417]. 

[20] N. Arkani-Hamed, Y. Grossman, and M. Schmaltz, Split fermions in extra dimensions and 
exponentially small cross-sections at future colliders, Phys. Rev. D61 (2000) 115004, 
[hep-ph/9909411]. 

[21] S. B. Giddings and S. D. Thomas, High energy colliders as black hole factories: The end of 
short distance physics, Phys. Rev. D65 (2002) 056010, [hep-ph/0106219]. 

[22] S. Dimopoulos and G. L. Landsberg, Black holes at the LHC, Phys. Rev. Lett. 87 (2001) 
161602, [hep-ph/0106295]. 

[23] P. C. Argyres, S. Dimopoulos, and J. March-Russell, Black holes and sub-millimeter 
dimensions, Phys. Lett. B441 (1998) 96-104, [hep-th/9808138]. 



- 18 - 



[24] D. M. Eardley and S. B. Giddings, Classical black hole production in high-energy collisions, 
Phys. Rev. D66 (2002) 044011, [gr-qc/0201034]. 

[25] U. Sperhake, V. Cardoso, F. Pretorius, E. Berti, and J. A. Gonzalez, The high-energy 
collision of two black holes, Phys. Rev. Lett. 101 (2008) 161101, [arXiv: 0806 . 1738]. 

[26] M. Shibata, H. Okawa, and T. Yamamoto, High-velocity collision of two black holes, Phys. 
Rev. D78 (2008) 101501, [arXiv: 0810. 4735]. 

[27] U. Sperhake et al., Cross section, final spin and zoom-whirl behavior in high- energy black 
hole collisions, Phys. Rev. Lett. 103 (2009) 131102, [arXiv: 0907. 1252]. 

[28] M. W. Choptuik and F. Pretorius, Ultra Relativistic Particle Collisions, arXiv: 0908. 1780. 

[29] C. M. Harris, P. Richardson, and B. R. Webber, CHARYBDIS: A black hole event generator, 
JEEP 08 (2003) 033, [hep-ph/0307305]. 

[30] M. Cavaglia, R. Godang, L. Cremaldi, and D. Summers, Catfish: A Monte Carlo simulator 
for black holes at the LEG, Comput. Phys. Commun. 177 (2007) 506-517, [hep-ph/0609001]. 

[31] D.-C. Dai et al, BlackMax: A black-hole event generator with rotation, recoil, split branes 
and brane tension, Phys. Rev. D77 (2008) 076007, [arXiv : 0711 . 3012]. 

[32] J. A. Frost et al, Phenomenology of Production and Decay of Spinning Extra- Dimensional 
Black Eoles at Eadron Colliders, JEEP 10 (2009) 014, [arXiv: 0904. 0979]. 

[33] R. C. Myers and M. J. Perry, Black holes in higher dimensional space-times, Ann. Phys. 172 
(1986) 304. 

[34] W. G. Unruh, Second quantization in the Kerr metric, Phys. Rev. D10 (1974) 3194-3205. 

[35] G. W. Gibbons, Vacuum polarization and the spontaneous loss of charge by black holes, 
Commun. Math. Phys. 44 (1975) 245-264. 

[36] P. Candelas, P. Chrzanowski, and K. W. Howard, Quantization of electromagnetic and 
gravitational perturbations of a Kerr black hole, Phys. Rev. D24 (1981) 297-304. 

[37] A. C. Ottewill and E. Winstanley, The renormalized stress tensor in Kerr space-time: 
General results, Phys. Rev. D62 (2000) 084018, [gr-qc/0004022]. 

[38] M. Casals and A. C. Ottewill, Canonical quantization of the electromagnetic field on the Kerr 
background, Phys. Rev. D71 (2005) 124016, [gr-qc/0501005]. 

[39] D. Ida, K.-y. Oda, and S. C. Park, Rotating black holes at future colliders: Greybody factors 
for brane fields, Phys. Rev. D67 (2003) 064025, [hep-th/0212108]. 

[40] C. M. Harris and P. Kanti, Eawking radiation from a (4+n)-dimensional black hole: Exact 
results for the Schwarzschild phase, JEEP 10 (2003) 014, [hep-ph/0309054]. 

[41] C. M. Harris and P. Kanti, Eawking radiation from a (4+n)- dimensional rotating black hole, 
Phys. Lett. B633 (2006) 106-110, [hep-th/0503010]. 

[42] D. Ida, K.-y. Oda, and S. C. Park, Rotating black holes at future colliders. II: Anisotropic 
scalar field emission, Phys. Rev. D71 (2005) 124039, [hep-th/0503052]. 

[43] G. Duffy, C. Harris, P. Kanti, and E. Winstanley, Brane decay of a (4 +n)- dimensional 
rotating black hole: Spin-0 particles, JEEP 09 (2005) 049, [hep-th/0507274]. 

[44] M. Casals, P. Kanti, and E. Winstanley, Brane decay of a (4+n)- dimensional rotating black 
hole. II: Spm-1 particles, JEEP 02 (2006) 051, [hep-th/0511163]. 



- 19 - 



[45] V. Cardoso, M. Cavaglia, and L. Gualtieri, Black hole particle emission in higher- dimensional 
spacetimes, Phys. Rev. Lett. 96 (2006) 071301, [hep-th/0512002]. 

[46] V. Cardoso, M. Cavaglia, and L. Gualtieri, Hawking emission of gravitons in higher 
dimensions: Non- rotating black holes, J 'HEP 02 (2006) 021, [hep-th/0512116]. 

[47] D. Ida, K.-y. Oda, and S. C. Park, Rotating black holes at future colliders. Ill: Determination 
of black hole evolution, Phys. Rev. D73 (2006) 124022, [hep-th/0602188]. 

[48] M. Casals, S. R. Dolan, P. Kanti, and E. Winstanley, Brane decay of a (4+n)-dimensional 
rotating black hole. Ill: Spin- 1/2 particles, JHEP 03 (2007) 019, [hep-th/0608193]. 

[49] M. Casals, S. R. Dolan, P. Kanti, and E. Winstanley, Bulk emission of scalars by a rotating 
black hole, JHEP 06 (2008) 071, [arXiv : 0801 . 4910]. 

[50] J. M. Bardeen, W. H. Press, and S. A. Tcukolsky, Rotating black holes: Locally nonrotating 
frames, energy extraction, and scalar synchrotron radiation, Astrophys. J. 178 (1972) 347. 

[51] S. Chandrasekhar, The mathematical theory of black holes, . Oxford, UK: Clarendon (1992) 
646 p. 

[52] A. L. Dudley and J. D. Finley, Separation of Wave Equations for Perturbations of General 
Type-D Space-Times, Phys. Rev. Lett. 38 (1977) 1505-1508. 

[53] K. D. Kokkotas, Quasinormal modes of the Kerr-Newman black hole, Nuovo Cim. B108 
(1993) 991-998. 

[54] E. Berti and K. D. Kokkotas, Quasinormal modes of Kerr-Newman black holes: Coupling of 
electromagnetic and gravitational perturbations, Phys. Rev. D71 (2005) 124008, 
[gr-qc/0502065]. 

[55] M. Rogatko and A. Szyplowska, Massive Fermion Emission from Higher Dimensional Black 
Holes, Phys. Rev. D79 (2009) 104005, [arXiv: 0904. 4544]. 



- 20 - 



